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Abstract. After a short historical review, we present four popular substructuring methods: 
FETI-1, BDD, FETI-DP, BDDC, and derive the primal versions to the two FETI methods, called P- 
FETI-1 and P-FETI-DP, as proposed by Fragakis and Papadrakakis. The formulation of the BDDC 
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1. Introduction. Substructuring methods are among the most popular and 
widely used methods for the solution of systems of linear algebraic equations obtained 
by finite element discretization of second order elliptic problems. This paper provides 
a review of recent results on the equivalence of several substructuring methods in a 
common framework, complemented by some details not published previously. 

We first give a brief review of the history of these methods (Section [2]). After 
introducing the basic concepts of substructuring (Section [3]), we formulate the dual 
methods, FETI-1 and FETI-DP (Section and derive their primal versions, r- 
FETI-1 and P-FETI-DP, originally introduced in [20]. However the derivation was 
omitted in [20]. Next, we formulate the primal methods, BDD and BDDC (Section[5]). 
Finally, we study connections between the methods in Section [6l We revisit our recent 
proof that the P-FETI-DP is in fact the same method as the BDDC [35] and the 
preconditioner by Cros ^Sj. Next, we translate some of the abstract ideas from [TOll^U] 
into a framework usual in the domain decomposition literature. We recall from [20] 
that for a certain variant of FETI-1, the P-FETI-1 method is the same algorithm as 
BDD. Then we derive a recent abstract result by Fragakis [19] in this special case 
to show that the eigenvalues of BDD and that particular version of FETI-1 are the 
same. It is notable that this is the variant of FETI-1 devised to deal with difficult, 
heterogeneous problems [jLJ . 

2. Historical remarks. In this section, we provide a short overview of iterative 
substructuring, also known as non-overlapping domain decomposition. Rather than 
attempting a complete unbiased survey, our review centers on works connected to the 
BDD and FETI theory by the second author and collaborators. 

Consider a second order, selfadjoint, positive definite elliptic problem, such as 
the Laplace equation or linearized elasticity, discretized by finite elements with 
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characteristic element size h. Given sufficient boundary conditions, the global stiffness 
matrix is nonsingular, and its condition number grows as O {h~^) for ft, — > 0. However, 
if the domain is divided into substructures consisting of disjoint unions of elements 
and the interior degrees of freedom of each substructure are eliminated, the resulting 
matrix on the boundary degrees of freedom has a condition number that grows 
only as O [H~^h~^) , where H >> h is the characteristic size of the substructure. 
This fact has been known early on (Keyes and Gropp [22]); for a recent rigorous 
treatment, see Brenner [3]. The elimination of the interior degrees of freedom is 
also called static condensation, and the resulting reduced matrix is called the Schur 
complement. Because of the significant decrease of the condition number, one can 
substantially accelerate iterative methods by investing some work up front in the 
Choleski decomposition of the stiffness matrix on the interior degrees of freedom 
and then just run back substitution in each iteration. The finite element matrix 
is assembled separately in each substructure. This process is called subassembly. 
The elimination of the interior degrees of freedom in each substructure can be done 
independently, which is important for parallel computing: each substructure can be 
assigned to an independent processor. The substructures are then treated as large 
elements, with the Schur complements playing the role of the local stiffness matrices 
of the substructures. See [22, 43^ for more details. 

The process just described is the background of primal iterative sub structuring 
methods. Here, the condition that the values of degrees of freedom common to several 
substructures coincide is enforced strongly, by using a single variable to represent 
them. The improvement of the condition number from O to 0[H~^h~^^, 

straightforward implementation, and the potential for parallel computing explain 
the early popularity of iterative substructuring methods [22j . However, further 
preconditioning is needed. Perhaps the most basic preconditioner for the reduced 
problem is a diagonal one. Preconditioning of a matrix by its diagonal helps to take 
out the dependence on scaling and variation of coefhcients and grid sizes. But the 
diagonal of the Schur complement is expensive to obtain. It is usually better to avoid 
computing the Schur complement explicitly and only use multiplication by the reduced 
substructure matrices, which can be implemented by solving a Dirichlet problem on 
each substructure. Probing methods (Chan and Mathew [6]) use such matrix- vector 
multiplication to estimate the diagonal entries of the Schur complement. 

In dual iterative substructuring methods, also called FETI methods, the condition 
that the values of degrees of freedom common to several substructures coincide 
is enforced weakly, by Lagrange multipliers. The original degrees of freedom are 
then eliminated, resulting in a system for the Lagrange multipliers, with the system 
operator consisting essentially of an assembly of the inverses of the Schur complements. 
Multiplication by the inverses of the Schur complements can be implemented by 
solving a Neumann problem on each substructure. The assembly process is modified 
to ensure that the Neumann problems are consistent, giving rise to a natural coarse 
problem. The system for the Lagrange multipliers is solved again iteratively. This 
is the essence of the FETI method by Farhat and Roux [TB], later called FETI-1. 
The condition number of the FETI-1 method with diagonal preconditioning grows 
as O (h^^) and is bounded independently of the number of substructures (Farhat, 
Mandel, and Roux [T7|). For a small number of substructures, the distribution of the 
eigenvalues of the iteration operator is clustered at zero, resulting in supcrconvergence 
of conjugate gradients; however, for more than a handful of substructures, the 
supcrconvergence is lost and the speed of convergence is as predicted by the O (h^^) 
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growth of the condition number [17j . 

For large problems and large number of substructures, asymptotically optimal 
preconditioners are needed. These preconditioners result typically in condition 
number bounds of the form O (log" (1 + H/h)) (the number 1 is there only to avoid the 
value log 1 = 0). In particular, the condition number is bounded independently of the 
number of substructures and the bounds grow only slowly with the substructure size. 
Such preconditioners require a coarse problem, and local preconditioning that inverts 
approximately (but well enough) the diagonal submatrices associated with segments 
of the interfaces between the subtructures or the substructure matrices themselves. 
The role of the local preconditioning is to slow down the growth of the condition 
number as ft. ^ 0, while the role of the coarse problem is to provide global exchange 
of information in order to bound the condition number independently of the number 
of substructures. Many such asymptotically optimal primal methods were designed 
in the 1980s and 1990s, e.g., Bramble, Pasciak, and Schatz [2113, Dryja [11], Dryja, 
Smith, and Widlund T3], Dryja and Widlund T4], and Widlund [46]. However, 
those algorithms require additional assumptions and information that may not be 
readily available from finite element software, such as an explicit assumption that 
the substructures form a coarse triangulation and that one can build coarse linear 
functions from its vertices. 

Practitioners desire methods that work algebraically with arbitrary substructures, 
even if a theory may be available only in special cases (first results on extending 
the theory to quite arbitrary substructures are given in Dohrmann, Klawonn, and 
Widlund [10] and Klawonn, Rheinbach, and Widlund [23]). They also prefer methods 
formulated in terms of the substructure matrices only, with minimal additional 
information. In addition, the methods should be robust with respect to various 
irregularities of the problem. Two such methods have emerged in early 1990s: the 
Finite Element Tearing and Interconnecting (FETI) method by Farhat and Roux [18] , 
and the Balancing Domain Decomposition (BDD) by Mandel [31]. Essentially, the 
FETI method (with the Dirichlet preconditioner) preconditions the assembly of the 
inverses of the Schur complements by an assembly of the Schur complements, and the 
BDD method preconditions assembly of Schur complements by an assembly of the 
inverses, with a suitable coarse problem added. Of course, the assembly weights and 
other details play an essential role. 

The BDD method added a coarse problem to the local Neumann-Neumann 
preconditioner by DeRoeck and Le Tallec [[41 , which consisted of the assembly (with 
weights) of pseudoinverses of the local matrices of the substructure. Assembling the 
inverses of the local matrices is an idea similar to the Element-by-Element (EBE) 
method by Hughes et al. [21]. The method was called Neumann-Neumann because 
the preconditioner requires solution of Neumann problems on all substructures, in 
contrast to an earlier Neumann-Dirichlet method, which, for a problem with two 
substructures, required the solution of a Neumann problem on one and a Dirichlet 
problem on the other [46]. The coarse problem in BDD was constructed from the 
natural nuUspace of the problem (constant for the Laplace equation, rigid body 
motions for elasticity) and solving the coarse problem guaranteed consistency of local 
problems in the preconditioner. The coarse correction was then imposed variationally, 
just as the coarse correction in multigrid methods. The O (log^ (1 -I- H/Kf) bound was 
then proved |31| . 

In the FETI method, solving the local problems on the substructures to eliminate 
the original degrees of freedom has likewise required working in the complement of 
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the nuUspace of the substructure matrices, which gave a rise to a natural coarse 
problem. Since the operator employs inverse of the Schur complement (solving a 
Neumann problem) an optimal preconditioner employs multiplication by the Schur 
complement (solving a Dirichlet problem), hence the preconditioner was called the 
Dirichlet preconditioner. The O (log^ (1 + H/Kj) bound was proved by Mandel and 
Tezaur and O (log^ (1 + H/h)) for a certain variant of the method by Tezaur fli] . 
See also Klawonn and Widlund [25^ for further discussion. 

Because the interface to the BDD and FETI method required only the 
multiplication by the substructure Schur complements, solving systems with the 
substructure Schur complements, and information about the substructure nullspace, 
the methods got quite popular and widely used. In Cowsar, Mandel, and Wheeler [7], 
the multiplications were implemented as solution of mixed problems on substructures. 
However, neither the BDD nor the FETI method worked well for 4th order problems 
(plate bending). The reason was essentially that both methods involve "tearing" 
a vector of degrees of freedom reduced to the interface, and, for 4th order problems, 
the "torn" function has energy that grows as negative power of h, unlike for 2nd order 
problems, where the energy grows only as a positive power of logl//i. The solution 
was to prevent the "tearing" by fixing the function at the substructure corners; then 
only its derivative along the interface gets "torn" , which has energy again only of the 
order \og\/h. Preventing such "tearing" can be generally accomplished by increasing 
the coarse space, since the method runs in the complement to the coarse space. For 
the BDD method, this was relatively straightforward, because the algebra of the 
BDD method allows arbitrary enlargement of the coarse space. The coarse space 
that does the trick contains additional functions with spikes at corners, defined by 
fixing the value at the corner and minimizing the energy. With this improvement, 
O {\og^ {1 + H/h)) condit ion number bound was proved and fast convergence was 
recovered for 4th order problems (Le Tallec, Mandel, and Vidrascu [28l [29]). In 
the FETI method, unfortunately, the algebra requires that the coarse space is made 
of exactly the nullspace of the substructure matrices, so a simple enlargement of 
the coarse space is not possible. Therefore, a version of FETI, called FETI-2, 
was developed by Mandel, Tezaur, and Farhat [3H], with a second correction by 
coarse functions concentrated at corners, wrapped around the original FETI method 
variationally much like BDD, and the O (log'^ (1 + H/h)) bound was proved again. 
However, the BDD and FETI methods with the modifications for 4th order problems 
were rather unwieldy (especially FETI-2), and, consequently, not as widely used. 

The breakthrough came with the Finite Element Tearing and Interconnecting - 
Dual, Primal (FETI-DP) method by Farhat et al. [TS], which enforced the continuity 
of the degrees of freedom on a substructure corner as in the primal method by 
representing them by one common variable, while the remaining continuity conditions 
between the substructures are enforced by Lagrange multipliers. The primal variables 
are again eliminated and the iterations run on the Lagrange multipliers. The 
elimination process can be organized as solution of sparse system and it gives rise to 
a natural coarse problem, associated with substructure corners. In 2D, the FETI-DP 
method was proved to have condition number bounded as O (log^ (1 -I- H/Kf) both for 
2nd order and 4th order problems by Mandel and Tezaur [37]. However, the method 
does not converge as well in 3D and averages over edges or faces of substructures 
need to be added as coarse variables for fast convergence (Klawonn, Widlund, and 
Dryja Farhat, Lesoinne, and Pierson [IS]), and the O (log^ {1 + H/h)) bound 
can then be proved t27J. 
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The Balancing Domain Decomposition by Constraints (BDDC) was developed 
by Dohrmann [9j as a primal alternative the FETI-DP method. The BDDC method 
uses imposes the equality of coarse degrees of freedom on corners and of averages 
by constraints. In the case of only corner constraints, the coarse basis functions are 
the same as in the BDD method for 4th order problems from [28l [29] . The bound 
O (log^ (1 + Hlh)) for BDDC was first proved by Mandel and Dohrmann ^33]. The 
BDDC and the FETI-DP are currently the most advanced versions of the BDD and 
FETI families of methods. 

The convergence properties of the BDDC and FETI-DP methods were quite 
similar, yet it came as a surprise when Mandel, Dohrmann, and Tezaur |34j proved 
that the spectra of their preconditioned operators are in fact identical, once all the 
components are same. This result came at the end of a long chain of ties discovered 
between BDD and FETI type method. Algebraic relations between FETI and BDD 
methods were pointed out by Rixen et al. [40], Klawonn and Widlund [25], and 
Fragakis and Papadrakakis [20j . An important common bound on the condition 
number of both the FETI and the BDD method in terms of a single inequality 
was given Klawonn and Widlund [25]. Fragakis and Papadrakakis [20], who derived 
certain primal versions of FETI and FETI-DP preconditioners (called P-FETI-1 and 
P-FETI-DP), have also observed that the eigenvalues of BDD and a certain version 
of FETI are identical along with the proof that the primal version of this particular 
FETI algorithm gives a method same as BDD. The proof of equality of eigenvalues of 
BDD and FETI was given just recently in more abstract framework by Fragakis [l9] . 
Mandel, Dohrmann, and Tezaur [33] have proved that the eigenvalues of BDDC and 
FETI-DP are identical and they have obtained a simplified and fully algebraic version 
(i.e., with no undetermined constants) of a common condition number estimate for 
BDDC and FETTDP, similar to the estimate by Klawonn and Widlund [25] for BDD 
and FETI. Simpler proofs of the equality of eigenvalues of BDDC and FETI-DP were 
obtained by Li and Widlund [30| . and by Brenner and Sung 5,, who also gave an 
example when BDDC has an eigenvalue equal to one but FETI-DP does not. A primal 
variant of P-FETI-DP was proposed by Cros [8', giving a conjecture that P-FETI-DP 
and BDDC is in fact the same method, which was first shown on a somehow more 
abstract level in our recent work 35J. 

It is interesting to note that the choice of assembly weights in the BDD 
preconditioner was known at the very start from the work of DeRoeck and Le Tallec 
[4T| and before, while the choice of weights for FETI type method is much more 
complicated. A correct choice of weights is essential for the robustness of the methods 
with respect to scaling the matrix in each substructure by an arbitrary positive 
number (the "independence of the bounds on jumps in coefficients"). For the BDD 
method, such convergence bounds were proved by Mandel and Brezina [32J, using a 
similar argument as in Sarkis [42^ for Schwarz methods; see also Dryja, Sarkis, and 
Widlund [12]. For the FETI methods, a proper choice of weights was discovered only 
much later - see Rixen and Farhat [39], Farhat, Lesoinne and Pierson [16] for a special 
cases, Klawonn and Widlund [25| for a more general case and convergence bounds, 
and a detailed discussion in Mandel, Dohrmann, and Tezaur [34] . 

3. Substructuring Components for a Model Problem. We first show how 
the spaces and operators we will work with arise in the standard substructuring theory 
for a model problem obtained by a discretization of the second order elliptic problem. 
Consider a bounded domain C decomposed into nonoverlapping subdomains 
(alternatively called substructures) denoted fi^, i = 1, A^, which form a conforming 
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triangulation of the domain Q. Each substructure is a union of a uniformly bounded 
number of Lagrangean PI or Ql finite elements, such that the nodes of the finite 
elements between substructures coincide. The boundary of Qi is denoted by dfli. The 
nodes contained in the intersection of at least two substructures are called boundary 
nodes. The union of all boundary nodes of all substructures is called the interface T 
and Ti is the interface of substructure fi^. The space of vectors of local degrees of 
freedom on Ti is denoted by Wi and W = Wi x ■ ■ ■ x Wn- Let Si : Wi —>■ Wi be the 
Schur complement operator obtained by eliminating all interior degrees of freedom of 
ili, i.e., those that do not belong to interface Ti. We assume that the matrices Si are 
symmetric positive semidefinite and consider global vectors and matrices in the block 
form 



Wi 



WN 



weW, S 



'Si 



Sn 



(3.1) 



The problem we wish to solve is the constrained minimization of energy, 
— a {u, u) — (r, u) — > min subject to u G W, 



(3.2) 



where C is the space of all vectors of degrees of freedom on the substructures 
that coincide on the interfaces, and the bilinear form 

a (u, v) — {Su, v) , Vu, w € W, 

is assumed to be positive definite on W. In the variational form, problem (13. 2|) can 
be written as 

ueW : a{u, v) = (r, v) , e W. (3.3) 
The global Schur complement S :W ^ W' associated with a is defined by 

a{u,v) ^ (Su,v) , yu,veW. (3.4) 



Defining R as the natural embedding of the space W into the space W, i.e. 



R:W ^W, R:ueW< 



ueW, 



(3.5) 



we can write (|3.3p equivalently as the system of linear algebraic equations 

Su = r, where S = R^SR. (3.6) 

The BDDC and FETI-DP, as the two- level preconditioners, are characterized by 
the selection of certain coarse degrees of freedom, such as values at the corners and 
averages over edges or faces of substructures (for their general definition see, e.g., [26]). 
So, we define W C W a.s the subspace of all functions such that the values of any 
coarse degrees of freedom have a common value for all relevant substructures and 
vanish on dfl, and such that 



W CW CW. 
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The space W has to be selected in the design of the preconditioner so that the bihnear 
form a (■, ■) is positive definite on W. The operator S : W ^ W associated with a is 
defined by 



a{u, v) = (^Su, , Vu, V G W. 



Remark 3.1. The idea to restrict the bilinear form a(-, •) from the space W into 
the subspace W is closely related to the concept of subassembly as employed in iSOf . 
In formulation of dual methods from the FETI family, we introduce the matrix 

B ~ [Bi, . . . , Bn] , 

which enforces the continuity across substructure interfaces and it is defined as follows: 
each row B corresponds to a degree of freedom common to a pair of substructures i 
and j. The entries of the row are zero except for one +1 in the block i and one —1 in 
the block j, so that the condition 

Bu = ueW, 

and using (|3.5[) . clearly 

BR = 0. (3.7) 

An irnportant ingredient of substructuring methods is the averaging operator E : 
W defined as 

E^R^Dp, (3.8) 

where Dp : W W is a given weight matrix such that the decomposition of unity 
property holds, 

ER = I. (3.9) 

In terms of substructuring, E is an averaging operator that maps the substructure 
local degrees of freedom to global degrees of freedom. 

The last ingredient is the matrix Bp) constructed from B as 

Bd ~ [DpiBi, . . . , DpinBn] , 



where the matrices Dpi are determined from Dp, see [271 [34] for details. 
Finally, we shall assume, cf., e.g., [31| eq. (10)], that 

BJ)B + RE = I, (3.10) 

which easily imphes EBj,B = E{I - RE) = E - ERE = 0, and so 

B'^BdE^ = Q. (3.11) 

4. P-FETI Family of Methods. We review the FETI-1 and FETI-DP 
preconditioners followed in each case by a formulation of their primal versions denoted 
as P-FETI- 1 and P-FETI-DP, respectively 
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4.1. P-FETI-1. In the case of the FETI-1 method, the problem (I32D is 
formulated as minimization of total subdomain energy subject to the continuity 
condition 

ifl {w, w) — (/, w) min subject to w &W, Bw = 0, (4.1) 

which is equivalent to a saddle point system: find [w^X) G x A such that 

Sw + S^A = /, 

Bw ^ 0. ^ ' 

First, note that S is invertible on nuUi? and A is unique up to a component in nuUi?"^, 
so A is selected to be range -B. Let Z be matrix with Unearly independent columns, 
such that 

range Z = nuUS*. (4.3) 
Since S is semi-definite, it must hold for the first equation to be solvable that 

/ - B^\ e ranges' = (null 5)^ = (range Z)^ = nuUZ'^, 
so, equivalently, we require that 

Z'^if - B'^X) = 0. (4.4) 

Eliminating w from the first equation of (|4.2p as 

w^S+{f-B'^X) + Za, (4.5) 

substituting in the second equation of (|4.2p and rewriting (|4.4p . we get 

BS+B^X - BZa = BS+f, 
^Z^B^X = -Z^f. 

Denoting G — BZ and F — BS^B^ this system becomes 

FA - Ga = BS+f, 
-G^X = -Z^f. ^^-^^ 

Multiplying the first equation by (Cf^QG) ^ G^'^Q, where Q is some symmetric and 
positive definite scaling matrix, we can compute a as 

a = (G^QG)"^ G^Q{FX - BS+f). (4.7) 

The first equation in (|4.6|) thus becomes 

FX-G {G^QGy^ G^Q{FX ~ BS+f) = BS+f. (4.8) 

Introducing 

P = I-QG{G'^QG)'^G^, (4.9) 



as the Q-orthogonal projection onto nullG"'", we get that (|4.8p corresponds to the 
first equation in (|4.6p muhiphed by . So, the system (I4.6P can be written in the 
decoupled form as 

P^FA = P^BS+f, 
G^X = Z^f. 

The initial value of A is chosen to satisfy the second equation in (14. 6p , so 

Ao = QGiG^QCy^Z^f. (4.10) 

Substituting Aq into (|4.7p gives initial value of a as 

ao = (G^QG) G^Q(FAo - BS+f). (4.11) 

bmce we are looking for A £ null G^, the FETI-1 method is a preconditioned conjugate 
gradient method applied to the system 



P^ FPX = P^BS+f. 



(4.12) 



In the primal version of the FETI-1 preconditioner, the assembled and averaged 
solution u is obtained from (|4.5|) . using equations (|4.11|) and (|4.10p . as 



u — Ew 
= E[S+{f ^B^\^) + Zao\ 

S+{f - B^Xo) + Z {G^QGy' G^QiFXo - BS+f) 

S+{f - S^Ao) + Z (G'^QGy^ G^Q{BS+B^Xo - BS+f) 

I-Z {G^QGy^ G^Qb] S+{f - B^Xo) 

J-Z {G^QGy^ G'^QB^ S+ (/ - B^QG {G'^QGy^ Z^ 



= E 
^ E 
= E 
= E 



E'r 



= EH'^S+HE'^r 

= Mp^FETir, 

where we have denoted by 



H = I - B^QG{G'^QG) ^ Z^ , 



and so 



Mp_FETI = EH^S+HE'^, 
is the associated primal preconditioner P-FETI-1, same as J2(T, eq. (79)]. 



(4.13) 



(4.14) 



4.2. P-FETI-DP. In the case of the FETI-DP, the problem ([321) is formulated 
as minimization of total subdomain energy subject to the continuity condition 



1 



a {w, w) — (/, w) min subject to w G W, Bw — 0. 



(4.15) 



Compared to the formulation of FETI-1 in (|4.ip . we have now used the subspace 
W C W such that the operator S associated with a (•, •) on the space W is positive- 
definite. In this case, (14.15^ is equivalent to setting up a saddle point system: find 
{w,X) &W X A such that 

Sw + B^X = f, , . 

Bw =0. ^ ' ' 

Since S is invertible on W, solving for w from the first and substituting into the 
second equation of (|4.16p . we get 

BS-^B^X = BS-^f, (4.17) 

which is the dual system to be solved by preconditioned conjugate gradients, with the 
Dirichlet preconditioner defined by 

MpETi-DP = BdSB]j. (4.18) 

Next, we will derive the P-FETI-DP preconditioner using the original paper by 
Farhat et. al. [H] in order to verify the P-FETI-DP algorithm given in gOl eq. (90)] 
for the corner constraints. We split the global vector of degrees of freedom u into the 
vector of global coarse degrees of freedom denoted by Uc and the vector of remaining 
degrees of freedom denoted by Ur- We note that we could perform a change of basis, 
cf., e.g., [miinilSn] to make all primal constraint (such as averages over edges or faces) 
explicit, i.e., each coarse degrees of freedom would correspond to an explicit degree of 
freedom in the vector Uc- Thus, we decompose the space W as, cf. [Ml Remark 5], 

W = Wc®Wr, (4.19) 

where the space Wc consists of functions that are continuous across interfaces, have 
a nonzero value at one coarse degreejaf freedom at a time and zero at other coarse 
degrees of freedom, and the space Wr consists of functions with coarse degrees of 
freedom equal to zero^ The solution splits into the solution of the global coarse 
problem in the space Wc and the solution of independent subdomain problems on the 
space Wr- 

Let i?^*^ be a map of global coarse variables to its subdomain component, i.e.. 



( R^' \ 



let Br be an operator enforcing the interface continuity of by 

BrUr =0, Br={ bP . . . S^^^ 



and let the mappings E'^ and distribute the primal residual r to the subdomain 
forces and to the global coarse problem right-hand side, respectively. 

The equations of equilibrium can now be written, cf. [151 eq. (9)-(10)], as 

SrrWr 

2^ He Ore Wr + 2^ Uc Occ -Kc Wc — Jc, 

2=1 1=1 

E Bi'^wi'^ = 0, 



=1 
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where the first equation corresponds to independent subdomain problems, second 
corresponds to the global coarse problem and the third enforces the continuity of 
local problems. This system can be re- written as 



{SrcRc) 
Br 



Sc. 








(4.20) 



where fr = E^r, fc — E^r, and the blocks are defined as 



N 










/ c.(1)r(1) 
Ore Ji-c 


\ 


Ucc — / , ^"-c '-'cc^'-c J 








1 Src^c — 










[ 


Srr j 




\ Ore Ji-c 





Remark 4.1. Note that the system \4.2(]^ is just the expanded system ^.16^ . 
Expressing Ur from the first equation in (|4.20p . we get 

Ur = S^r^ [fr — SrcRcUc — BJ\) . 

Substituting for Ur into the second equation in (|4.20p gives 

Sqc^c ~ {SrcRc) Srr ^ ^ fc ^ {SrcRc) Srr fr: 

where 5**^, — Sec — R^ Sj^Srr SrcRc- Inverting S*^, we get that 



5** 



fc {SrcRc) Srr fr + (SrcRc) Srr ^r ^ 



After initialization with A = 0, which [5D1 [TO] does not say, but it can be used, cf., 
e.g., [ini Section 6.4], the assembled and averaged solution is 

U — ErUr + EcUc 

= ErSrr ^fr SrcRcS cc (^fc ~ {SrcRc) Srr /'"^ ^ ^" 
+ EcS*e ^/c — (SrcRc) Srr fr 

ErSrr fr ErSrr SrcRcS cc fc 

+ ErSrr SrcRcS cc (SrcRc) Srr fr^ 

+ EcS^e fc — EcS^f, (SrcRc) Srr fr 
= ErSrr fr + 

+ (Ec — ErSrr SrcRc) S ^.^ (^fc ~ (SrcRc) Srr fr 
— ^Ip-FETI-Dpr, 



where 



Mp-FETI-DP — ErSrr Er + 



(4.21) 



^" (^Ec ErSrr SrcRc^ Sec ("^c Rc SrcSrr Er ) 

is the associated preconditioner P-FETI-DP, same as [lOl eq. (90)]. 
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5. BDD Family of Methods. We recall two primal preconditioners from 
the Balancing Domain Decomposition (BDD) family by Mandel in [31]; namely 
the original BDD and Balancing Domain Decomposition by Constraints (BDDC) 
introduced by Dohrmann [9 . 

5.1. BDD. The BDD is a Neumann-Neumann algorithm, cf., e.g., [2], with 
a simple coarse grid correction, introduced by Mandel The name of the 

preconditioner comes from an idea to balance the residual. We say that f G is 
balanced if 

Z^E^v = 0. 

Let us denote the "balancing" operator as 

C = EZ, (5.1) 

so the columns of C are equal to the weighted sum of traces of the subdomain zero 
energy modes. Next, let us denote by ScS the S — orthogonal projection onto the 
range of C, so that 

and by Pc the complementary projection to ScS, defined as 

Pc^I- ScS. (5.2) 
The BDD preconditioner |311 Lemma 3.1], can be written in our settings as 



M. 



BDD 



s- 



-1 



(/ - Sc'S^ ES+E'^S{I - ScS) + ScS 

J - Sc'S^ ES+E'^{SS-^ - SScSS-'^) + ScSS-^ 
= PcES+E^P^ + Sc (5.3) 

where Sc serves as the coarse grid correction. See [31] (32] , and [5D] for details. 

5.2. BDDC. Following a similar path as Li and Widlund [30], we will assume 
that each constraint can be represented by an explicit degree of freedom and that we 
can decompose the space W as in (|4.19p . We note that the original BDDC in [^ [33] 
is mathematically equivalent, but algorithrnically it treats the corner coarse degrees 
of freedom and edge in the definition of W in different ways. The BDDC is the 
method of preconditioned conjugate gradients for the assembled system (j3.6l) with 
the preconditioner Mbddc defined by cf. [3D1 eq. (27)], 

MbdDC — Tsub + Tq, 

where Tgub = Ej.S~^Ej is the subdomain correction obtained by solving independent 
problems on subdomains, and To = E'i! {^'^ S'^) 'i'^E^ is the coarse grid 
correction. Here ^ are the coarse basis functions defined by energy minimization, 

min . 



Since we assume that each constraint corresponds to an explicit degree of freedom, 
the coarse basis functions ^' can be easily determined via the analogy to the discrete 
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harmonic functions, discussed, e.g., in [321 Section 4.4]; are equal to 1 in the coarse 
degrees of freedom and have energy minimal extension with respect to the remaining 
degrees of freedom Ur, so they are precisely given as 



^ = 

Then, we can compute 



Rc 



\T,T q\T, _ ( f>T _nTcTc-l ^ ( ^cc Sj^ \ f Rc 



followed by 



Sec S^^Sj,^ S-rc^C 



— {^Eq EfS^.^ Si^qRc^ ^ cc {"^c Rq S rr ) * 

So, the BDDC preconditioner takes the form 

Mbddc = ErS-^Ej+ (5.4) 

+ {Ec - ErS:^^SrcRc) s*c {eJ ~ sj^s:^^^ eJ) . 

6. Connections of the Preconditioners. We review from [20l Section 8] that 
a certain version of P-FETI-1 gives exactly the same algorithm as BDD. Next, we state 
the equivalence of P-FETI-DP and BDDC preconditioners. Finally, we translate the 
abstract proof relating the spectra of primal and dual preconditioners [191 Theorem 
4] in the case of FETI-1 and BDD. 

Theorem 6.1 ( [20. Section 8]). If Q is chosen to be the Dirichlet preconditioner, 
the P-FETI-1 and the BDD preconditioners are the same. 

Proof. We wiU show that the P-FETI-1 in (|4T4|) with Q = BdSB^ is the same 
as the BDD in So, similarly as in [20, pp. 3819-3820], from (|4T3l) we get 

H^I- B^QG {G^QCy^ 

= I~ B'^BdSB'^BZ {z'^b'^BdSb'^bzY^ Z^ 

= I-An{Z^A„)-'z^, 

where 

A/j = B^ BjjSB^BZ. 
Using ((3T01) . definitions of C in jsll), S in dSJ]), and because S'Z by gj]), 

Ar = (/ - E'^R^) S{I - RE) Z 

^ SZ - SREZ - E'^R^SZ + E'^R^SREZ 
^SZ- SRC - E'^R^SZ + E'^SC 

= {e'^s- spdj c, 
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and similarly 



Z'^Ar = Z^ [e^S - SRj C 
= C^SC - Z^SREZ 
= C^SC. 

Using the two previous results, (|5.2p and symmetries of S and Sc, we get 



HE' = I-Ar Z'Aj,) Z' )E 



E^-Ar(Z^Ar) 'z^e^ 



^E' ^ E'S-SR]C C'SC] C 



^E'^ - (^E'^S - S, 



'c 



^E^ - E^SSc + SRSc 
= E'^ [l- SSc) + SRSc 
= E'^P^ + SRSc. 
Next, the matrix Sc satisfies the relation 

ScR^ SS^ SRSc — ScR^ SRSc — ScSSc 

c(c^5c) ^c'^sc(c^sc^ ^C^ 

= C(C^5C) ^C^ = Sc. 

Because by definition PcG = 0, using p.9p we get for some Y that 

PcES+SRSc = PcE [I + ZY) RSc 

= PcERSc + PcEZYRSc 

5c - 5c - 0, 

and the same is true for the transpose, so Sc R^ S S^ E'^ P^ = 0. 

Using these results, the P-FETI-1 preconditioncr from (|4.14l) becomes 

A'lp—pETi ^ EH^ S^HE^ 

= (ScR^S + PcE) S+ {e'^P^ + SRSc) 
= ScR^SS+E^P^ + ScR^SS+SRSc 

+ PcES+E^P^ + PcES+SRSc 
= PcES+E^P^ + Sc, 

14 



and we see that (|6.ip is the same as the definition of BDD in (|5.3p . □ 

Theorem 6.2. The P-FETI-DP and the BDDC preconditioners are the same. 

Proof. The claim follows directly comparing the definitions of both 
preconditioners, P-FETI-DP in eq. (|4?2T|) and the BDDC in eq. □ 

Corollary 6.3. Comparing the preconditioner proposed by Cros 18, eq. 4-S] with 
the definitions {4.21^ and \5.4% it follows that this preconditioner can be interpreted 
as either, P-FETI-DP or BDDC. 

In the remaining, we will show the equality of eigenvalues of BDD and FETI-1, 
with Q being the Dirichlet preconditioner. 

Lemma 6.4. The two preconditioned operators can be written as 



MfetiT^ {BdSBD [bS+B'^ 
MbddS = (eS+eA (R^SR) , 



where 

S+ = H^S+H. 

Proof. First, MpETi = B^SB]^, which is the Dirichlet preconditioner. From 
(|4.12p . using the definition of H by (|4 13[) . we get 

= p^bs+b'^p 

= {I- G{G'^QG)-^G^Q'^) BS+B'^ (/ - QG{G'^QG)-'^G'^) 
= {B- BZ{G^QG)-^G^Q^B) 5+ {B^ - B^QG{G^QGy^ B^) 
= B{I- Z{G'^QG)-^G'^QB) S+ (/ - B^QG{G^QG)-^Z^) B^ 
= BH^S+HB^ = BS+B'^. 

Next, S is defined by (|3.6p . By Theorem 16. 1[ we can use (|4.14p for Mbdd to get 
Mbdd = EH^S+HE^ = ES+E^. □ 

Before proceeding to the main result, we need to prove two technical Lemmas 
relating the operators S and S~^. The first Lemma establishes [lH Assumptions (13) 
and (22)] as well as ^ Lemma 3] for FETI-1 and BDD. 

Lemma 6.5. The operators S, defined by liS.l]) and Theorem \6.4\ resp., satisfy 

S+SR = R, (6.2) 

S+SS+ = S+. (6.3) 
Moreover, the following relations are valid 

BS+SR = 0, (6.4) 

S+B^BdSS+E^ = 0. (6.5) 
Proof. First, from (|4.3p and symmetry of S it follows that 

HS=(^I- B^QG (G'^QCy^ Z^) 5* 
= S- B'^QG {G'^QCy^ Z'^S = S. 
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Using =I-Z {C^QG) ^ C^QB we get 

H'^S+S ^H^ {I + ZY) ^H'^ + H^ZY 

I - Z [G^QCy^ G'^QB 



ZY 



+ ZY - Z (G^QG) ^G^QGY 



= H'^ + ZY - ZY = iJ^, 



5+5* = H^S+HS = H^S+S = H'^ . 
Finally, from previous and (|3.7p . we get (|6.2p as 

S+SR = H^R = (l-Z {G^QGY^ G^Qb) R = R, 



and since H is a projection, we immediately get also 

5^+55+ = H^S+ = H^H^S+H = S+. 

Next, ([Q]) follows directly from noting ([SJ]) . 

Using (in21)-(IS31) and (iXTU)l - (PTT|l . we get (ES]) as 

S+B^BdSS+E^ = S+ {I- E^R^) SS+E^ 

= S+SS+E^ - S+E^R^SS+E^ 

= S+E^ - S+E'^R^E'^ 

= S+ {I - E^R^) E^ 

= S+B^BdE^ = 0. □ 

Next Lemma is a particular version of [El Theorem 4] for FETI-1 and BDD. 
Lemma 6.6. The following identities are valid: 



Td (Mfeti^) - [MbddS) Td, Td = ES+B' , 

Tp (^MbddS^ - {MfetiJ") Tp, Tp = (Mfeti^) BdSR. 

Proof. Using the transpose of (|6.5p and (|6.4p . we derive the first identity as 

Td {Mfeti^) = ES+B^BdSB^BS+B^ 

= ES+ {I - E'^R^) S{I- RE) S+B^ 

= £^5^+5 (/ - RE) S+B^ - ES+E^R^SS+B^ 

+ ES+E^R^SRES+B^ 
= ES+SBlBS+B^ - ES+E^R^SS+B^+ 
'eS+eA (R^SR) Td 



= i MbddS] Td. 
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Similarly, using (|6.5p and (|6.4p . we derive the second identity as 

Tp [MbddS^ ^ [Mfeti^F] BdSRES+ E'^ SR 

= {MfetiT) BdS {I - BlB) S+ {I ~ B^Bd) SR 
= {Meet I J") BdSS+ (I - B^Bd) SR 

- {MfetiT) BdSB]jBS+SR 

+ {MfetiT) BdSBIbS+B^BdSR 
= MfetiBS+ B^ BdSS+ R^ SR 

- [MfetiT) BdSB]jBS+SR 

+ {MfetiT) (BdSBI) {bS+B'^^ BdSR 

= {MfetiJ") {Mfeti^F) BdSR. 
= {Mfeti:F) Tp. □ 

Theorem 6.7. Under the assumption of Lemma \6.6[ the spectra of the 
preconditioned operators MbddS and Mfeti-i^^ satisfy the relation 

o (MbddS) \ {1} = a {Mfeti-iT) \ {0, 1} . 

Moreover, the multiplicity of any common eigenvalue A ^ 0, 1 is identical for the two 
preconditioned operators. 

Proof. Let up) he a (nonzero) eigenvector of the preconditioned FETI-1 operator 
corresponding to the eigenvalue Xp. Then, by Lemma 16.61 we have 

Tp {Mfeti-i^)up = (^AIpppS^ Tpup, 

so Tpup is an eigenvector of the preconditioned BDD operator corresponding to the 
eigenvalue \p^ provided that Tpup ^ 0. So assume that Tpup = 0. But then it is 
also true that 

= BpSR(Tpup) = BpSRES+B^up 
= BpS (/ - B^B) S+B'^up = BpSS+B^up - BpSBlBS+ B'^up 
= BpSS+B'^up - [Mfeti^F) up = BpSS+B^up - Xpup, 

so 

BpSS'^ B^up = Xpup. 
Note that, by dO]) and dST]), we get 

(BdSS+B'^'^^ = BpSS+B^BpSS+B'^ 

= BpSS+ (/ - SS+B'^ 
= BpSS+SS+B'^ - BpSS+E'^R^SS+B^ 
= BpSS+B^ - BpSS+E^R^B'^ 
= BpSS+B'^, 
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so BdSS^ is a projection and therefore — 0, 1. 

Next, Let up be a (nonzero) eigenvector of the preconditioned BDD operator 
corresponding to the eigenvalue Ap. Then, by Lemma 16.61 we have 

Tp [MbddS) = {MFETir)Tp, 

so Tpup is an eigenvector of the preconditioned FETLl operator corresponding to 
the eigenvalue Ap, provided that Tpup ^ 0. So assume that Tpup = 0. But then 
also using (|6.2p and (|3.9p . we get 

{) = Td [Tpup) = Td {MfetiT) BdSRup 

- (MbddS) TdBdSRup - (^MbddS) ES+B^BdSRup 

= MbddSES+ (/ - S'i?up 

= MbddSES+SRup - MbddSES+ E^ R^ SRup 

= MbddSup - MbddSES+ E^ R^ SRup 

= MbddSup - [aIbddS^ up, 
which is the same as 

Xpup — X^pUp = Ap (1 — Ap) Up ~ 0, 

and therefore Ap = 0, 1. ^ 

Finally, let A 7^ 0, 1 be an eigenvalue of the operator MbddS with the multiplicity 
m. From the previous arguments, the eigenspace corresponding to A is mapped by the 
operator Tp into an eigenspace of Mpeti-i^ and since this mapping is one-to-one, 
the multiplicity of A corresponding to Mpeti-i^ is n > m. By the same argument, 
we can prove the opposite inequality and the conclusion follows. □ 
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